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Abstract We present predictions for the evolution of FTEs generated by localized bursts of reconnection 
on a planar magnetopause that separates a magnetosheath region of high densities and weak magnetic 
field from a magnetospheric region of low densities and strong magnetic field. The magnetic fields present 
a shear angle of 1 05 degrees. Reconnection forms a pair of FTEs each crossing the magnetopause in the field 
reversal region and bulging into the magnetosphere and magnetosheath. At their initial stage they can be 
characterized as flux tubes since the newly reconnected magnetic field lines are not twisted. Reconnection 
launches Alfvenic perturbations that propagate along the FTEs generating high-speed jets, which move the 
pair of FTEs in opposite directions. As the FTE moves, it displaces the ambient magnetic field and plasma 
producing bipolar magnetic field and plasma velocity signatures normal to the nominal magnetopause in 
the regions surrounding the FTE. The combination of the ambient plasma with the FTE flows generates a 
vortical velocity pattern around the reconnected field lines. During its evolution the FTE evolves to a flux 
rope configuration due to the twist of the magnetic field lines. The alfvenic perturbations propagate faster 
along the part of the FTE bulging into the magnetosphere than in the magnetosheath, and due to the differ- 
ences between the plasma and magnetic field properties the perturbations have slightly different signatures 
in the two regions. As a consequence, the FTEs have different signatures depending on whether the satellite 
encounters the part bulging into the magnetosphere or into the magnetosheath. 


1. Introduction 

Flux transfer events (FTEs) have been interpreted as the result of transient and localized reconnection occur- 
ring at the dayside magnetopause [Russell and Elphic, 1978]. They occur frequently [Rijnbeek et al., 1984] and 
are believed to be responsible for a significant fraction of the flux transport from the dayside to nightside 
magnetosphere [ Cowley , 1984]. In satellite data, they are characterized by a positive-negative (standard FTE) 
or a negative-positive (reverse FTE) pulse in the component of the magnetic field perpendicular to the nom- 
inal magnetopause ( B N ). This pulse is thought to result from the draping of the ambient magnetic field lines 
around the reconnected flux tube combined with the formation of an axial field-aligned current along the 
flux tube itself [Paschmann et al., 1 982]. The field-aligned current along the FTE is associated with a twist in 
the field and allows the interpretation of FTEs as a flux rope instead of a flux tube [Paschmann eta!., 1982]. 
The total field magnitude inside the FTE can be characterized by an overall peak or a decrease (crater), the 
latter being suggested as an initial stage of FTE evolution [Zhang et al., 2010]. The decrease in the total field 
in crater FTEs consists of two "trenches" in the magnetic field at the edges of the crater, with a slight increase 
in the field in the center reflecting the existence of a weak force-free helical current similar to that observed 
in flux ropes [LaBelle et al., 1987]. Support for the hypothesis that FTEs are flux ropes instead of flux tubes is 
also given in Zhang et al. [2008] where magnetic field observations from multiple satellites were fitted to 
a flux rope model. Other features appearing to be invariant characteristics of FTE are a mixture of mag- 
netospheric and magnetosheath particle populations and an increase in the total pressure inside the FTE 
compared to that outside the event [Paschmann et al., 1 982], 

Once FTEs form, they move across the magnetopause plane, away from the reconnection site, in response 
to pressure gradient and curvature forces [Cooling et al., 2001]. Dailey etal. [1985] investigated a series of 
23 FTE events and determined an average £ x 8 drift velocity of 1 25 km/s for the reconnected flux tube. 
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In investigating the spatial properties of FTEs using ISEE 1 and 2 satellite data, Saunders etal. [1984] reported 
the first evidence for plasma vorticity in or around FTEs tubes. They associated this vorticity with magnetic 
field twisting and a torsional Alfven wave propagating antiparallel to the ambient field. Liu et al. [2008] 
found, using Time Fiistory of Events and Macroscale Interactions during Substorms (TFIEMIS) multiple 
satellite data, that flow perturbations outside of FTEs in the magnetosphere are approximately dipolar. 

Korotova et al. [2009] developed a flow signature taxonomy for FTEs. For events moving faster than the 
ambient media the flow perturbation directly upstream and downstream from the events should be in the 
direction of the event motion. Flows on the flanks should be opposite to event motion. They should also be 
accompanied by bipolar inward/outward (outward/inward) flow perturbations normal to the nominal mag- 
netopause in the magnetosphere (magnetosheath). Events moving with the ambient flow should generate 
no flow perturbations in the ambient plasma. Alfven waves propagating parallel (antiparallel) to the FTE's 
axial magnetic field should generate anticorrelated (correlated) flow perturbations within the core region 
of FTEs. 

Farrugia etal. [1987] developed a simple incompressible flow model for the interaction of the ambient field 
and fluid with a rigid flux tube of cylindrical geometry. The flow obtained was a classical potential flow 
(V x v = 0), and with a knowledge of the flow field from hydrodynamics, the magnetic field could also 
be determined. Their results show that plasma approaching the tube is decelerated and deflected at the 
periphery of the tube. The flow has a dipolar form in the rest frame of the flux tube, and the magnetic field 
perturbations obtained from this flow resemble FTE signatures. The normal component executes a bipolar 
excursion about zero, while the tangential components suffer substantial changes which maximize when B N 
passes through zero. 

Recently, Zhang et al. [201 1] investigated the flow properties around 3701 FTEs observed between May and 
October in both 2007 and 2008 by TFIEMIS spacecraft. The flow perturbations observed in conjunction with 
the magnetic signatures show that FTEs moving in an antisunward direction on the magnetopause are asso- 
ciated with flow vortices just inside the magnetopause. In order to understand the flow perturbations they 
have developed a 2-D MFID model that starts from the solution obtained by Farrugia et al. [1 987] but remove 
the assumption that the flow is inviscid and impose finite viscosity near the FTE. They conclude that viscosity 
plays a key role in the formation of flow vortices in the low-latitude boundary layer, and when it is consid- 
ered the simulation reproduces the observations qualitatively. However, they do not discuss the mechanism 
responsible for the increase in viscosity close to the FTE nor investigate the evolution of magnetic field and 
plasma inside the FTE. 

Ku and Sibeck [1997, 1998a, 1998b, 2000] performed a series of 2-D magnetohydrodynamic (MHD) simula- 
tions to investigate the signatures of FTEs produced by a Petschek-type reconnection occurring close to 
the subsolar point. In their first work [Ku and Sibeck, 1 997] reconnection is started by a step function time- 
dependent localized increase in resistivity. The plasma is at rest at the beginning of the simulation, and the 
initial magnetosheath and magnetosphere field and plasma parameters are chosen in a way to reproduce 
realistic ratios observed by satellites for the magnetosheath-magnetosphere system. They observe that a 
bulge of merged field lines develops almost entirely on the magnetosheath side of the magnetopause 
and grows with time. The bulge moves along the magnetopause at a fraction of the magnetosheath 
Alfven speed. A plasma jet in the direction of the bulge's motion forms in a boundary layer on the magne- 
tospheric side of the magnetopause. This jet is on magnetic field lines with depressed magnetic field 
strengths and approach and slightly exceeds the magnetosheath Alfven velocities. Outside the bulge a gen- 
tler reverse flow appears on the magnetosheath side. This reverse flow is associated with a vortical structure 
that follows the bulge during its motion. Ku and Sibeck [1998a] found that events moving opposite to a 
background magnetosheath flow slow down but intensify, whereas events moving in the direction of the 
magnetosheath flow accelerate but weaken. The vortical flow outside the bulge on the magnetosheath is 
weakened for events moving in the direction of the magnetosheath flow, whereas it is strengthened for 
events moving opposite to the flow. They associate the FTEs signature with the vortical flow at the magne- 
tosheath. Ku and Sibeck [1 998b] checked for the effects of bursty reconnection and secondary reconnection 
by considering a background resistivity and turning on and off the enhanced local resistivity. In this case 
after the enhanced resistivity is switched off, secondary reconnection sites appear on the trailing edge of 
the event where thin current sheets are formed. Merging at these sites inflates and propels the main FTE. 
They also fill a magnetic island with hot tenuous plasma near the subsolar point. Transient magnetic islands 
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also form on the trailing edge of the main FTE, forming secondary FTEs, but they remain small and later 
disappear. Besides giving important information about FTE evolution and the associated perturbation 
in the ambient plasma and magnetic field, the 2-D simulations performed in Ku and Sibeck [1 997, 1 998a, 

1 998b, 2000] only consider perfectly antiparallel fields, a particular case that is rarely observed at the mag- 
netopause. Also, in their 2-D simulations the FTE evolution is restricted to the plane perpendicular to the 
magnetopause as well as the evolution of associated perturbations. The combination of these two facts 
leads to the interpretation of the FTE as the bulge forming close to the field reversal region and caused 
difficulties in separating magnetospheric from magnetosheath evolution. 

Global simulations to investigate the generation mechanisms of FTEs have been performed by Fedder et al. 
[2002], Raeder [2006], and Dorelli and Bhattacharjee [2009]. Fedder et al. [2002] reported that the FTE is 
formed via global separator reconnection. They invoke a current-driven instability as the driver of recon- 
nection, and as separator reconnection proceeds, strong parallel currents form along the reconnected flux 
tube, resulting in a twist of the field lines. Raeder [2006] had investigated the generation mechanism of FTEs 
during strongly southward interplanetary magnetic field. In his simulations FTEs are generated by multiple 
X line reconnection. The multiple X lines form sequentially for periods of significant dipole tilt. This gives a 
seasonally modulated formation rate for FTEs, with no FTE formation during equinox and formation in the 
winter hemisphere during solstices. The model explains the 8 min quasiperiodicity by the convection and 
reformation time of flux ropes at the magnetopause. In the Dorelli and Bhattacharjee [2009] simulation as 
the subsolar magnetopause current density grows in time, the local geometry of the magnetic field changes 
in such a way that the stagnation point of the flow moves away from the subsolar point, producing a large 
flow shear at the subsolar magnetopause. The subsolar flow shear is unstable to the formation of flow vor- 
tices, which grow and coalesce as the FTE grows at the subsolar magnetopause. The flow vortices within 
the FTE drive 3-D magnetic reconnection resulting in a complex magnetic topology. In performing global 
simulations, the authors focused mainly on the generation mechanisms of FTEs and not on the 3-D evolu- 
tion of reconnected magnetic field lines and associated plasma and magnetic field perturbation. Also, as 
stressed by Dorelli and Bhattacharjee [2009], reconnection in the case of Fedder et al. [2002] and Raeder 
[2006] are due numerical resistivity, and the results obtained in the simulation by Raeder [2006] and Dorelli 
and Bhattacharjee [2009] for plasma temperature and density inside the FTEs are not in agreement with 


In this work we extended the investigations performed in Ku and Sibeck [1 997, 1 998a, 1 998b, 2000] by using 
a local 3-D resistive MFID simulations and considering a sheared magnetic field configuration. Similar work 
has been performed by Otto [1 990, 1 995], Ma et al. [1 995], La Belle-Fiamer et al. [1 995], and Otto et al. [1 995]. 
Flere we focus on the 3-D nature of the flow in and around the FTE during its evolution and the differences 
in the FTE signature and related perturbation in the magnetosheath and magnetosphere plasma and mag- 
netic field. In section 2 we present the model equations, initial, and boundary conditions. In section 3 we 
present the results obtained in our numerical experiment. Finally, in section 4, we summarize our findings 
and draw our conclusions. 

2. Model Setup 

The model used to perform this work was successfully applied to study the evolution of plasma and mag- 
netic field in the solar atmosphere [ Santos and Buchner, 2007; Santos et al., 2008, 2011a, 201 1 b]. It solves an 
appropriate set of resistive MFID equations (in normalized units) 
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Table 1. Values Used for the Normalization of the Macroscopic 
Variables 


Macroscopic Variable 

Normalization Value 

B 

B 0 = 1 x 10“ 8 T 

L 

L 0 =4x 10 5 m 

T 

T 0 = 10 6 K 

P 

p o = Po£- o =3.6x10-' n Pa 

P 

=2.2x10-21 kg 

zk b' 0 

V 

v a = v A = _ = 191 , 642 m/s 

u \/wo 

t 

f 0 = ^ = 2.1 s 
°0 


together with Ohm's law, Ampere's law, and taking into account the ideal gas equation for a fully ionized 
plasma 


E = -ux B + rjj 

(5) 

V x 8 = nj 

(6) 

p = 2riK B T 

(7) 


Here p is the plasma mass density, u is the plasma velocity, 8 is the magnetic field, p is the plasma pressure, 

T is the plasma temperature, p 0 is the normalization value used for the plasma beta, and y is the polytropic 
index (y = /). Table 1 presents the values used in the normalization of the macroscopic variables. 

The system of equations is solved in an equidistant cartesian grid (131 x 131 x 131 grid points) that covers 
a 3-D volume of 1 2 x 24 x 32 Mm 3 . This gives a grid resolution of approximately 91 x 183 x 244 km in x, y 
and, z directions, respectively. The equations are solved using a second-order scheme in time (Leapfrog) and 
a second-order scheme in space (centered differences). The simulation volume presents six boundaries: four 
lateral, top, and bottom boundaries. We apply Neumann boundary conditions on all boundaries with j = 0, 
apart on the lateral boundaries where we have considered the boundary closed for perpendicular flows. In 
this way we allow a free outflow on the top and bottom boundaries and a reflective flow on the other four 
boundaries. We keep V ■ 8 close to zero by applying the projection method for divergence cleaning [Toth, 
2000 ], 

The initial conditions for the plasma and magnetic field are the same used in Otto [1 990] and represent a pla- 
nar current sheet (the magnetopause) which separates regions of different plasma properties. These initial 
conditions are reproduced below: 


p{x) = p, + p 2 tanh 



( 8 ) 


p x (x) = p 00 + 2b z0 A fa tanhx + (1 - /? 0 )A 2 cosh 2 x 


(9) 


b y (x) = yjb 2 ya + /? 0 A 2 cosh 2 x 


( 10 ) 


b z (x) = b z0 - A b tanhx 


( 11 ) 
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Figure 1 . Time evolution of the magnetic field (tubes) and plasma flow (arrows) for simulation of FTE formation via local- 
ized Petschek reconnection, (top to bottom) The snapshots show results at f = 210 s, 294 s, and 378 s. The view is in the 
positive x direction and macroscopic variables are presented in normalized units. 
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The free parameters 6, b 5b , and p 00 determine 
the asymptotic magnetic shear, the amplitude 
of the magnetic field in the magnetosheath 
(normalized units), and the typical thermal 
pressure (normalized units). The additional 
parameters p v p 2 ,x 0 , and d x determine the 
jump in mass density (normalized units), the 
position, and width of the transition region 
where the jump occurs. Magnetosheath plasma 
is initially at rest since we consider reconnec- 
tion near the stagnation point. However, as 
stressed by Korotova et at. [2009], the flow 
results obtained in this situation can also be 
applied to the situation when the FTE moves 
faster than the ambient plasma. As in Otto 
[1 990], the y and z directions do not correspond 
to the y and z directions described in GSE or 
GSM coordinate systems. 

In Petschek's model [Petschek, 1964] dissipation 
occur within a small region inside the current 
sheet, known as the diffusion region. In the dif- 
fusion region, disruption of the current leads 
to the breakdown of the frozen-in approxima- 
tion allowing reconnection to occur. In view of 
the fact that the physical process giving rise to 
reconnection, which is believed to be micro- 
scopic in nature, is still a matter of debate, we 
consider reconnection as initiated by a spot 
of localized increase in magnetic diffusivity 
[Buchner and Elkina, 2005, 2006]. As in Otto 
[1 990] we assume that the localized magnetic 
diffusivity is present for f > 0 and is given by 

r){x,y,z) = r] a cosh -2 - cosh -2 - cosh -2 - — — . 

3 ( 1 2 ) 

where r] a is the anomalous magnetic diffusivity. 


3. Simulation Results 

We perform a numerical experiment to inves- 
tigate the properties of FTEs formed by a 


Figure 2. Time evolution of the magnetic field 
(tubes) and flow field lines for simulation of FTE 
formation via localized Petschek reconnection, 
(top to bottom) The snapshots show results at 
f = 210 s, 294 s, and 378 s. The view is in the positive 
x direction. 
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localized Petschek-type reconnection in a 3-D space and its influence in the ambient plasma and magnetic 
field in the course of its evolution. This mechanism for generating FTEs goes back to the mechanism sug- 
gested by Russell and Elphic [1 978]. The merging occurs at a magnetopause marked by sharp transitions in 
density, temperature, and magnetic field strength and direction. The ratio of magnetosheath to magneto- 
spheric densities is 4, the ratio of field strengths is 0.75, and that of temperatures is 0.385. The magnetic field 
is subject to a shearing of approximately 1 05° close to the plane defined by x = 0 (magnetopause). The 
magnetosheath region is located atx > 0, while the magnetosphere is located atx < 0. 

The merging occurs due to a localized increase in the magnetic diffusivity (anomalous diffusivity, Buchner 
and Elkina [2005, 2006]) in a region around the center of the simulation box, known as diffusion region, 
where the normalized values of the diffusivity increase from zero to a peak value of 0.05 in normalized units. 
The increase in magnetic diffusivity causes a localized breakdown of the frozen-in approximation, allowing 
the field to change connectivity. This region also acts as a localized source of reconnection electric field and 
accompanying disturbances, as well as a generator of magnetohydrodynamic (MHD) waves. The MHD waves 
propagate rapidly outward from the dissipation region, along the reconnected field lines, and dominate 
energy conversion on large scale. 

Figure 1 presents a 3-D visualization of plasma flow (arrows) and magnetic field (tubes) patterns at three dif- 
ferent times following the onset of reconnection (f = 210 s, 294 s, and 378 s). The magnetic field lines start 
on points along aline going from z = 20 0 to z = 60L 0 , at the plane x = 0. The length of the arrows are pro- 
portional to the magnitude of the flow velocities. The color code depicts the velocity magnitude saturated 



-15 -10 -5 0 5 10 15 

Figure 3. Lateral view of magnetic field (tubes) and flow field lines at f = 378 s. The view is in the positive y direction. 
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Figure 4. Vortical flow around FTE 1 at f = 378 s. The magnetic field (tubes) is colored with the x component of the field, 
in normalized units, while the velocity vectors (arrows) are colored by the velocity magnitude. 

at 0.5v 0 , where v 0 is the normalization velocity (see Table 1), for better visualization. Reconnection takes 
place within a spherical region of radius « 5 L 0 , where L 0 is the normalization length scale (see Table 1 ), with 
center located at (x 0 ,y 0 ,z 0 ) = (0, 0,40) in normalized units. Two regions of merged field lines develop and 
move along the magnetopause, above, and below the diffusion region in the z direction. In the region above 
z = 40 plasma is accelerated mainly in the positive z direction, while in the region below z = 40 plasma 
is accelerated mainly in the negative z direction. These are the outflow jets characteristic of the reconnec- 
tion process. Figure la shows that at t = 210 s plasma is accelerated only in regions close to the diffusion 
region. Flowever, reconnection generates MHD waves which propagate rapidly outward from the dissipation 
region along the reconnected field lines. The velocity of propagation of these waves is larger in the magne- 
tosphere than in the magnetosheath, because the Alfven velocity is larger there. As a consequence, energy 
conversion is expected to take place faster along the part of the reconnected field lines that bulges into the 
magnetosphere than on the part that bulges into the magnetosheath. Figure 1 c shows that at f = 378 s 
the plasma jets have been formed along the part of the reconnected magnetic field lines that bulges into 
the magnetosphere, reaching the top and bottom of the simulation box. The jets do not form much further 
away from the diffusion region along the part of the reconnected field lines that bulges into the magne- 
tosheath. At t = 378 s the plasma velocity peaks at approximately 0.7v 0 , which corresponds to 1 34 km/s. 
This value is very close to the average value of 1 25 km/s found observationally by Dailey et at. [1 985] for the 
average £ x 6 drift velocity of the reconnected flux tubes. The maximum Alfven mach number obtained in 
the simulation volume at this time is about 1 .3 in a region close to the diffusion region. 



Figure 5. Magnetic field associated to the flux transfer event 1 (FTE 1) and flux transfer event 2 (FTE 2) at f = 378 s. 
Also shown is the density in a plane intersecting the two FTEs aty = 10t 0 . The color code depicts the density varia- 
tion in normalized units. 
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Associated with the plasma jets that form along the reconnected field lines, there is a vortical plasma flow. 
This vortical flow pattern is not visible in Figure 1 because it is very weak compared to the plasma jet gen- 
erated in the process. Figure3 shows the magnetic field and flow field lines at f = 210 s, 294 s, and 378 s. It 
is visible from this figure that the flow field lines gradually concentrate around the part of the reconnected 
field lines that bulges into the magnetosphere. Figure 3 shows a lateral view of Figure 3c. The vortical plasma 
flow is on the magnetospheric side, and it closes on the plasma jet. The vortical flow also appears around 
the part of the reconnected field lines that bulges into the magnetosheath but does not form much fur- 
ther away from the diffusion region. Figure 4 shows the velocity vectors in a plane at y = 10. We draw 
field lines starting in a region on that plane close to the center of the vortex and along a line in the x direc- 



Figure 6. (top) Component of the magnetic field and (bottom) plasma velocity perpendicular to the nominal magne- 
topause at a xz plane located at y = 10t 0 , for t = 378 s. The view is in the positive y direction, and macroscopic variables 
are presented in normalized units. 
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tion that crosses the vortex region in the plane y = 1 0. The field lines obtained are a combination of open 
(reconnected) and closed field lines. From the figure it seems that the strong jets are on closed field lines 
and the vortical structure is due the motion of the plasma around the reconnected field lines. This is simi- 
lar to the situation described in Korotova et al. [2009] and in the model of Farrugia et al. [1 987]. There is also 
a flow perpendicular to the diffusion region, along the x direction. This flow brings magnetic flux from the 
magnetosheath and magnetosphere to feed the reconnection process. 

We now analyze the spatial variation of the plasma parameters associated to the FTE. Different from the 2-D 
simulations performed by Ku and Sibeck [1997, 1998a, 1998b, 2000] or in cases where the magnetosheath 
and magnetospheric magnetic fields are perfectly antiparallel, here the FTE signature is not anymore asso- 
ciated to the bulge that forms around the field reversal region but it is distributed along the reconnected 
field lines as the Alfvenic perturbation propagates from the reconnection point. As a consequence, we can 
follow the evolution of the part of the FTE bulging into the magnetosphere and the part bulging into the 
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Figure 7. (top) Plasma density and (bottom) total magnetic field intensity at a xz plane located aty = 10L 0 , for f = 378 s. 
The view is in the positive y direction, and macroscopic variables are presented in normalized units. 
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magnetosheath separately and check for differences in FTE signatures. We select a region along a plane xz 
cutting the simulation volume at y = 10 (see Figure 5), for f = 378 s. This plane is approximately paral- 
lel to the direction of motion of the reconnected field lines and intersects the two regions of reconnected 
field lines formed above and below the diffusion region, in the z direction, which we call flux transfer event 
1 (FTE 1) and flux transfer event 2 (FTE 2), respectively. Figure 6 shows the components of the magnetic 
field (Figure 6a) and plasma velocity (Figure 6b) normal to the nominal magnetopause (x component). 
Figure 6a shows that both FTEs produce bipolar signatures in the normal component of the magnetic field. 
This bipolar signature is due to the motion of plasma around the FTE, which bends the ambient magnetic 
field lines, and due to the twist of the magnetic field inside the FTE. If FTE 1 passes a satellite located in a 
position atz > 0, the satellite will observe a positive-negative pulse in the normal component of the mag- 
netic field. This signature would correspond to a standard FTE in satellite data. If FTE 2 passes a satellite 
located in a position at z < 0, the satellite will see a negative-positive pulse in the normal component of 
the magnetic field. This signature would correspond to a reverse FTE in satellite data. Figure 6b shows that 
the FTEs are also associated with bipolar variations in the component of the velocity normal to the nomi- 
nal magnetopause. The velocity variations are anticorrelated to magnetic field variations and are caused by 
the plasma motion around the FTE. There is no normal component of the velocity associated with the FTE 
itself. Figure 6 also reveals that there is an asymmetry between the signatures of FTE 1 , which bulges into 
the magnetosphere, and FTE 2, which bulges into the magnetosheath, in the x components of velocity and 
magnetic field. This asymmetry is visible also in other plasma parameters. Figure 7 shows the plasma den- 
sity (Figure 7a) and magnetic field intensity (Figure 7b). The density in FTE 1 increases compared with the 
density in the ambient magnetosphere. This increase in density is associated with a decrease in the total 
magnetic field and an increase in temperature (not shown here). The density in FTE 2 decreases compared 
with the ambient magnetosheath density. This density decrease is associated with an increase in the total 
field and an increase in temperature (not shown here). The total field decreases in the trailing edge of FTE 
2, which is probably associated with the removal of magnetic flux from the reconnection region. The sig- 
nature in plasma density and magnetic field related to FTE 2 presents a kind of vortical structure, which is 
not visible around FTE 1 . From the results obtained here we conclude that the same FTE can have different 
signatures in satellite data, depending on whether the satellite encounters the part bulging into the mag- 
netosphere or into the magnetosheath. A plane xz cutting the simulation volume aty = -10 would show a 
mirror image of what we have described in Figures 6 and 7, with FTE 1 now penetrating the magnetosheath 
and FTE 2 the magnetosphere. Flowever, in this case, the bipolar variations in the velocity component nor- 
mal to the nominal magnetopause would be correlated to the magnetic field variations. Again, this velocity 
is caused by plasma motion around the FTE tube, and the normal component of the magnetic field is an 
effect of the bending of magnetic field lines around it. 

The distinct signatures from different trajectories relative to FTE 2 in the plane y = 1 0L 0 are presented in 
Figures 94 2. Figure 8 shows the four distinct trajectories T1 , T2, T3, and T4. These trajectories have a length 


Magnetosphere [ Magnetosheath 



T1 T2 T3 T4 


Figure 8. Plasma density near the magnetopause and four trajectories T1 , T2, T3, and T4 in the plane y = 1 0 perpen- 
dicular to the magnetopause. The perturbed magnetopause region corresponds to the presence of a FTE located on the 
magnetosheath side of the magnetopause. 


SANTOS ET AL. 


©2013. American Geophysical Union. All Rights Reserved. 


2019 



IgAGU Journal of Geophysical Research: Space Physics 


10.1 002/201 3JA01 8964 


of 1 6L 0 (« 6400 km) along the z direction, parallel to the magnetopause, going from z = 1 7L 0 to z = 33L 0 . 
Trajectory T1 is located at x = 1 L 0 , at 400 km from the center of the magnetopause into the magnetosheath. 
Trajectories T2, T3, and T4 are located at x = 2 L 0 (« 800 km), x - 4 L 0 (« 1 600 km), and x = 5 L 0 (a 2000 km), 
respectively. Figure 9 shows the variations along trajectory T4 of the component of the magnetic field 
normal to the magnetopause (Figure 9, top, left), the total field strength (Figure 9, top, right), and density 
(Figure 9, bottom), in normalized units. The component of the magnetic field normal to the magnetopause 
presents the typical bipolar signature characteristic of FTEs. Trajectory T4 starts at the magnetosheath and 
skims the FTE. When approaching the transition region between the magnetopause and magnetosheath, 
the total field increases, and the plasma density decreases. At the center of the FTE the total field strength 
reaches its maximum, and the plasma density reaches its minimum, returning to magnetosheath values 
at the end of the trajectory. Figurel 0 shows the same signatures for trajectory T3. This trajectory starts in 
the magnetosheath and crosses the transition region. When crossing the transition region, the total field 





Figure 9. Plots of the component of the magnetic field normal to (top, left) the magnetopause, (top, right) the total field 
strength, and (bottom) plasma density for trajectory T4. This trajectory have a length of 1 6Z_ 0 (« 6400 km) along the z 
direction, parallel to the magnetopause, going from z = 1 7L 0 to z = 33 L 0 , and is located at x = 5 L 0 (« 2000 km) from 
the center of the magnetopause into the magnetosheath. Thex axis corresponds to the length of the trajectory starting 
from z = 1 7L 0 . All plots are presented in normalized units. 



Figure 10. The same as Figure 9 but for trajectory T3 located at x = 4 L 0 (a 1600 km) from the center of the 
magnetopause into the magnetosheath. 
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strength (top, right) increases and the density (bottom) decreases, the peak values being larger than those 
obtained for trajectory T4. The results for trajectory T2 are presented in Figurel 1 . This trajectory starts at the 
magnetosheath and crosses the FTE itself. The total field strength (top, right) increases and reach its max- 
imum at the center of the FTE, where the component of the magnetic field normal to the magnetopause 
changes sign. At the trail part of the FTE the total field strength decreases to magnetosheath values, (bot- 
tom) The density decreases at the front part of the FTE, associated to the crossing of the magnetopause, and 
returns to increase inside the FTE itself reaching magnetosheath values. At the end of the trajectory the den- 
sity decreases again since it crosses the transition region. The total field strength and the component of the 
magnetic field normal to the magnetopause start to increase at the end of the trajectory when entering the 
transition region. Figure 12 shows the signatures for trajectory T1 .This trajectory starts at the magnetopause 
region and crosses the distorted magnetopause and transition region, returning to the magnetosheath at 
the end. The total field strength reaches its maximum at the center of the FTE and decreases below mag- 
netosheath values at the end of the trajectory. Plasma density monotonically increases when crossing the 
FTE and at the end of the trajectory, when entering the magnetosheath. It is important to notice that the 
trajectories are not perpendicular to the axis of the FTE. 



C 




Figure 1 1. The same as Figure 9 but for trajectory T2 located at x = 2 L 0 (« 800 km) from the center of the magnetopause 
into the magnetosheath. 





L 


Figure 12. The same as Figure 9 but for trajectory T1 located at x = 1E 0 (as 400 km) from the center of the magnetopause 
into the magnetosheath. 
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4. Summary and Conclusions 

Our results show that a localized Petschek-type reconnection generates two regions of reconnected field 
lines above and below the diffusion region in the z direction, which we have called flux transfer event 1 (FTE 
1) and flux transfer event 2 (FTE 2), respectively. At their initial stage the FTEs can be characterized as flux 
tubes since the newly reconnected magnetic field lines that form the FTE are not twisted. The FTEs move 
along the magnetopause plane, away from the reconnection site. During its motion the FTE displace the 
ambient plasma producing a bipolar signature in the component of the velocity normal to the nominal 
magnetopause. The plasma motion bends the ambient magnetic field producing a bipolar signature also in 
the component of the magnetic field normal to the nominal magnetopause. The bipolar magnetic field per- 
turbation is also present in the core of the FTE while there is no bipolar velocity signature visible there. This 
indicates that during its evolution the magnetic field lines that form the FTE are twisted, and the FTE evolves 
from a flux tube to a flux rope configuration. The magnetic field twist is associated with electric currents 
along the FTE. Strong plasma jets are observed in the core of the FTE, and the combination of the bipolar 
flow field around the FTE and the jets gives a vortical flow pattern around the FTE. 

Different from 2-D simulations or 3-D simulations where the magnetosheath and magnetospheric mag- 
netic fields are perfectly antiparallel, here the FTE signature is not anymore associated with the bulge 
that forms around the field reversal region but it is distributed along the reconnected field lines as 
the Alfvenic perturbation propagates from the reconnection point. As a consequence, we are able fol- 
low the evolution of the part of the FTE bulging into the magnetosphere and the part bulging into 
the magnetosheath separately and check for differences in FTE signature. We identify significant differ- 
ences between the signature of the FTE associated with the part that bulges into the magnetosphere 
and the part that bulges into the magnetosheath. In the part that bulges into the magnetosphere, 
the density increases inside the FTE compared to the ambient plasma. This density increase is asso- 
ciated with an increase in temperature, similar to the result obtained in the global simulations of 
Raeder [2006] and Dorelli and Bhattacharjee [2009]. The total magnetic field decreases inside the core of 
the FTE, characterizing a crater FTE. In the part that bulges into the magnetosheath, the density decreases 
inside the FTE compared to the ambient plasma. This density decrease is associated with an increase in 
temperature, as usually observed in FTEs. The total magnetic field also increases inside the FTE. 
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